Gaussian density fluctuations and Mode Coupling Theory for supercooled liquids 
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The equations of motion for the density modes of a fluid, derived from Newton's equations, 
are written as a linear generalized Langevin equation. The constraint imposed by the fluctuation- 
dissipation theorem is used to derive an exact form for the memory function. The resulting equations, 
solved under the assumption that the noise, and consequently density fluctuations, of the liquid are 
gaussian distrib uted, are equivalent to the random-phase-approximation for the static structure 
factor and to the well known ideal mode coupling theory (MCT) equations for the dynamics. This 
finding suggests that MCT is the canonical mean-field theory of the fluid dynamics. 



The description of the dynamics of supercooled liquids 
is one of the most intriguing goals in condensed-matter 
physics. Different approaches have been pursued in the 
last decades [0-0, originating both from the physics of 
liquids and from the physics of disordered systems. Many 
of these approaches strongly suggest that two different 
mechanisms for the decay of fluctuations are active in 
two different temperature ranges. A cross-over tempera- 
ture T c separates the region of 'weak' supercooling from 
the region of 'strong' supercooling. 

The description of the long-time dynamics and the 
associated evaluation of the transport coefficient in the 
strong supercooling region (below T c ) has proved an ex- 
tremely difficult task ^fli|]|]. The thermodynamic de- 
scription of the liquid state between T c and the calori- 
metric glass transition has been attempted, but no well 
defined connection between thermodynamics and dynam- 
ics has been achieved as yet. In the weak supercooling 
regime, detailed predictions for the space and time de- 
pendence of the long time decay of density correlations 
have been formulated using the ideal mode coupling the- 
ory (MCT) Q , one of the first approaches to identify the 
existence of the cross-over temperature. The agreement 
of MCT predictions with experimental findings H[l(| and 
exact correlation functions evaluated from molecular dy- 
namics simulations [pd|-|l3| — both for atomic and molec- 
ular models — supports the view that MCT is indeed 
able to describe the slow dynamics in weak supercooled 
states. 

Despite its remarkable practical success, the presence 
of apparently uncontrolled approximations in the deriva- 
tion of the MCT equations makes it difficult to gain in- 
sights into possible improvements of the theory. The aim 
of this Letter is to present a new derivation of the ideal 
MCT equations, starting from the microscopic equations 
for the evolution of the density (Newton's equations) and 
writing them as a linear generalized Langevin equation. 



A formally exact expression for the memory kernel is de- 
rived and, on making the approximation that the noise 
in the Langevin equation is Gaussian, the standard MCT 
equations are obtained. Note that the proposition of 
Gaussian noise implies that the density fluctuations are 
also gaussian [[l4| . The outcome of the present approach 
offers a route to MCT that involves clear and known 
approximations for static and dynamics. At the same 
time it suggests how MCT equations can be developed 
for more complex systems such as molecules and poly- 
mers. 

The density of a system composed of N particles lo- 
cated at positions r 3 is p(r,t) = YljLi — r jW)- The 
Fourier transform of p(r,t) is pk(t) = 2j=i e 4k ' rj In 
the case of a pairwise additive potential, the time evolu- 
tion of Pk(t) can be written as 

p k (t) = -^(k-r,-(t))V k ^«- 
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where is the Fourier transform of the pair potential. 
Eq.|l| is the Newton equation for the variables Pk(t). 

The equation of motion for the density (Eq. [[]) can be 
rewritten as, 

Pk{t) + to 2 kP k{t) = Tk{t) (2) 

where the linear term in P u_(t), that represents the ele- 
mentary excitations of the system, has been explicitly iso- 
lated. These excitations can be interpreted as 'phonons', 
that would oscillate with frequencies to\ if no interaction 
force amongst them were present. Combining Eqs.|l| and 
U the residual force ?\t{t) has the formal expression 

7v(t) = ulpk{t) - • r^)) V k ^ (t > - 
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The value of uif. that results in the least residual in- 
teraction between phonons, i.e. the one for which 

d(\^)/dui = o is mn 
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This choice for uj\ is also imposed by requiring the correct 
short-time limit of Eq.|l| p7| . The minimization principle 
produces an orthogonality condition between the force 
and the density variables at all times, 



(p-k(^k(t)> = 
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Substituting the expression for T k (t) (Eq. |) in this 
orthogonality condition leads to the exact Yvon-Born- 
Green (YBG) equation @] 
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where S^- = (/9- k/ o k )/./V is the static structure factor, 
/? = l/ksT, ks is the Boltzmann constant, T is the tem- 
perature and the symbols (■ ■ •} indicate equilibrium av- 
erages. 

Following the spirit of the Zwanzig-Mori (ZM) formal- 
ism fil| , the residual force can be written as sum of two 
contributions: (i) a term of dissipative origin, containing 
a memory function dependent on time. This represents 
slow deviations of the system from the integrable parts of 
phase space, (ii) A general random noise term that rep- 
resents a more 'chaotic' behavior for the system. Hence, 



^k(t) = - / 7k (* - t')d V (^{t')dt' + 
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where 7k (t) acts as memory function of the system and 
fk(t) as fluctuating force. Note that the previous choice 
for the residual force (Eq. ^) satisfies the orthogonality 
condition, Eq. ^[ 

The functions fk(t) and 7k (t) are not independent. In- 
deed, they must satisfy the fluctuation-dissipation theo- 
rem (FDT) to guarantee that the long time evolution 
of the system is consistent with the correct Boltzmann 
equilibrium distribution. This requires that the autocor- 
relation function of /k(i) is proportional to the memory 
function 7k (t). Also, the average over the noise is zero. 
Thus, 



</k(*)> = o 
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As a result of all these formal assignments the Newton's 
equations for the density fluctuations are rewritten in a 
form similar to a linear generalized Langevin equation as, 

0k(t) + w\ p k (t) + / 7k (t - t')dt<p k (t')dt' = /k(t) (10) 
Jo 

Note that the transition from the Newton's equations 
to the generalized Langevin equation is associated with a 
transition from averages over the initial conditions to av- 
erages over the realization of the noise (defined in Eqs. ||- 
|Io| ). FDT guarantees that the dynamics generated by 
this stochastic process leads to thermodynamic equilib- 
rium. 

Using Eq. [)| we calculate the explicit exact form for 
the memory function 20 as, 
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Eqs. |^-^T| constitute an exact definition of the noise 
process and have been obtained by using only FDT and 
the causality relation (p_ k (0)/ k (t)) =0. These equations 
have the same content as the equations derived in the ZM 
formalism p9| . To explicitly evaluate the 7k(i) function, 
approximations have to be made on the stochastic pro- 
cesses p-k{t) and /k(i), i-e. the density and the noise pj| . 

In this Letter, a simple approximation is considered. 
We assume that /k(^) is an additive gaussian process. 
This is the major approximation, since the assumption 
of a gaussian process and the linearity of Eq. 10 imply 
that Pk(t) is also a gaussian process. As shown in Eqs. || 
and [7] the true / k (t) is a non linear function of p k (i). This 
indicates that in an exact theory both quantities can not 
be simultaneously gaussian distributed Q . 

A less significant approximation consists in averaging 
over the velocities neglecting the fluctuations in the single 
particle kinetic energy, i.e. assuming 
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The gaussian nature of the noise ensures certain sim- 
plifications in the properties of pk{£) and this allows to 
calculate the multiple averages in Eq. Together with 
the approximation in Eq. [l2[ this allows to derive an 
expression for Jkit) which requires as input the density- 
density correlation functions. As a result, the gaussian 
approximation in the density gives, 
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The previous expression still contains the Fourier trans- 
form of the pair potential v k . As is well known |l7],|l8| , 
once the gaussian approximation for the density has 
been made, the YBG equation (Eq. ^J) can be consis- 
tently solved (2^] , providing the so-called Random-Phase- 
Approximation (RPA) |p7jp3| , c k = —(3v k , between the 
direct correlation function c k = (1 — l/S(k))/n and the 
potential v k . Eliminating in Eq.[13| v k in favor of c k one 
obtains, 
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7k (t) coincides with the memory function that has 
been originally calculated in Ref. [EJ within the ideal 
MCT framework. 

The reader may note that the choice of a gaussian pro- 
cess for p k (i), if used directly in Eq.||||, implies that the 
time evolution of the density-density correlation function 
is described by undamped harmonic modes [ p3|]2~|l with 
frequency uj k . In this approximation there is no inter- 
action between the modes and hence strictly speaking 
7k{t) = 0. In the present derivation, the memory kernel 
is not zero, since the assumption that the noise is gaus- 
sian is introduced only after the residual interactions are 
constrained to be a noise and dissipation and a formally 
correct expression for 7k(i) (Eq. ll]) is derived. This 
apparent contradiction is a consequence of the fact that 
the same approximations made before constraining phase 
space to have the Mori-type properties are not equiva- 
lent to those made after this fundamental constraint is 
forced upon the system. Therefore, the division of the 
residual forces between density waves into a dissipative 
and noise term is fundamental. This step, whilst not al- 
ways acceptable, is a conventional approach to dealing 



with complex many body forces, and the fact that the 
integrable motions have been clearly removed from the 
interactions prior to this division is quite satisfying. 

The present calculations show the following. 

(i) The assumption of gaussian properties for pk(t) and 
/k(i) allows to derive MCT. All approximations used in 
the conventional MCT are exact in this limit. Hence 
MCT can be seen, as a minimum, as an exact theory in 
the RPA limit and a fully consistent mean-field approx- 
imation to the dynamics of a complex system. Thus, 
the MCT dynamics is in the same class of universality 
as mean-field dynamics. This possibility was suggested 
some time ago on the basis of the analogies between the 
equations describing the schematic MCT models and the 
the dynamics of the order parameter in disordered p- 
spin models, solved under strict mean- field approxima- 
tion |||. 

(ii) The present approach offers a more direct inter- 
pretation to the uncontrolled approximation intrinsic in 
the conventional MCT. Indeed, as clearly stated in Ref. 
|jl| , a priori estimation of the quality of the conventional 
MCT are not known. The derivation of the basic equa- 
tions does not include any systematic expansion scheme. 
In the present approach, the more transparent approxi- 
mations suggest that improvement over the MCT predic- 
tions is feasible. In particular, it is possible to eliminate 
v k in Eq.[ll] by implementing higher order approxima- 
tions for the triplet correlation function. For example, 
it is possible to implement in the present scheme the 
Singwi-Sjolander (STLS) closure (^6|, which introduces 
controlled corrections to the gaussian statistics. Work is 
in progress on this topic. 

(iii) In the RPA, the equation for the memory kernel, 
Eq. [H], simplifies and the integral contributions can- 
cel out. This simplification does not occur in general. 
It may be shown that |27J these terms provide a better 
description of the short time dynamics as well as a renor- 
malization of the memory at infinite time. It can also be 
shown that these terms are dropped in the conventional 
MCT, when projecting onto the density pair subspace. 
In the present approach, these terms can be added to the 
conventional MCT result. This will provide a consistent 
expression for the short time memory function. 

(iv) An interesting aspect of the Gotze MCT deriva- 
tion is that it does not enforce the consistency condi- 
tions between static and dynamics and therefore the in- 
put structure factors do not have to be RPA in origin 
(and indeed, "exact" direct correlation functions calcu- 
lated from simulation data are often used as input in 
the theoretical calculations |ll|,Q). In the conventional 
MCT approach, v k disappears because of the gaussian 
approximation in calculating the normalization matrix 
of the projection [jlj. If consistency between statics and 
dynamics were enforced, the RPA would have followed, 
as discussed above. Moreover, Kawasaki M has recently 
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presented a derivation of MCT based on the (quadratic) 
density functional Ramakrishnan-Yussouf free energy of 
a liquid, where the effective interaction between the den- 
sity pairs is exactly — c(fc) / /3, suggesting that it is possible 
to correct the theory, without fundamentally affecting the 
universality class of its dynamics. Perhaps the answer to 
the intriguing question of why this is successful lies in 
formulating higher corrections to the YBG and memory 
kernel equations. 

(v) Finally, we stress that the present calculation will 
open new possibilities for an improved MCT description 
of the dynamics of deep supercooled liquids as well as 
extensions in order to deal with more complex systems, 
like molecules [ pilfjsl |30|] and polymers j3lj| . This is a 
strategy actively pursued, and it will be most interesting 
to see the outcome. 
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